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Abstract 

A complete study of the propagation of waves (namely longitudinal density and tem- 
perature waves, longitudinal and transversal velocity waves and heat waves) in turbulent 
superfluid helium is made in three situations: a rotating frame, a thermal counterflow, and 
the simultaneous combination of thermal counterflow and rotation. Our analysis aims to 
obtain as much as possible information on the tangle of quantized vortices from the wave 
speed and attenuation factor of these different waves, depending on their relative direction 
of propagation with respect to the rotation vector. 
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1 Introduction 

The most known phenomenological model, accounting for many of the properties of He II, 
given by Tisza pQ and Landau [2] is called the two-fluid model. The basic assumption is 
that the liquid behaves as a mixture of two fluids: the normal component with density p n and 
velocity v n , and the superfluid component with density p s and velocity v s . When the difference 

V := v n — v s between the normal and superfluid velocities, known also as counterflow velocity, 
exceed a certain critical velocity, a mutual friction F sn has to be included. This friction force 
is attributed to an interaction of the normal component with the vortices in the superfluid. 

Quantized superfluid vortices play an important role in the hydrodynamics of the fluid and 
they have been the object of many studies. The state of the fluid in which vortices are present, 
is referred to as the superfluid turbulent state. A review on superfluid turbulence can be found 
in Tough's paper [3] and in chapter 7 of Donnelly's book [3]. The quantized vortices created by 
applying a thermal counterflow form an irregular, spatially disordered tangle of lines. In this 
case, the vortex line density L (length of vortex line per unit volume) is Lh ~ J 2 V 2 , where 

V is the modulus of the relative velocity between the two components of the mixture and 7 
a temperature-dependent coefficient [3]. The vortex system is almost isotropic, provided that 
one neglects a small anisotropy induced by the imposed counterflow [5]. 
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The creation of the vortices cannot be made only in this way; in fact, the first studies of 
quantized vorticity involved a sample of He II rotating at constant angular velocity £1 exceeding 
a certain small critical value. The results brought to an ordered array of vortices aligned along 
the rotation axis, whose number density per unit area is given by Feynman's rule L/j = 2Q/k, 
where k = h/m = 9.97 10~ 4 cm 2 /sec is the quantum of circulation, with h Planck's constant 
and m mass of the helium atom. 

Now, an important question naturally arises: what happens if vortices are created by both 
rotation and counterflow? There has been only one experiment of which we are aware [6], on 
the formation of vortices in combined rotation and counterflow along the rotational axis. This 
experiment suggests that there exists a form of steady rotating turbulence, characterized by 
a vortex line density at given counterflow velocity V and angular velocity fi. Swanson et al. 
[6] found that at slow rotation the critical counterflow velocity above which the flow became 
turbulent was greatly reduced. The experimental observations showed that the two effects 
(thermal counterflow and rotation) are not merely additive, in fact for V high the measured 
values of L are always less than Lh + Lr- However, from our point of view the results of these 
experiments are purely qualitative because the authors didn't take the anisotropy of the vortex 
tangle in consideration which, as we will see in the last section of this paper, is essential to 
know the spatial distribution of the vortex tangle in liquid Helium II through measurements 
of second sound attenuation. 

The aim of this work is to study the propagation of longitudinal density and temperature 
waves, and longitudinal and transversal velocity waves and heat waves in the combined situa- 
tion of a rotating frame and of a cylindrical container in presence of thermal counterflow. The 
studies of the two separate cases of pure rotation and pure thermal counterflow are also con- 
sidered in order to give a more complete view of the wave propagation in these three different 
situations. The influence of the parameters characterizing the vorticity on the propagation of 
the waves is shown explicitly. The practical interest of this research is to obtain information 
on the vortex tangle from measurements on wave propagation. This is an important issue, 
because under the combined influence of rotation and counterflow the vortex tangle cannot be 
assumed isotropic. Then, we must find not only the vortex line density L but also the geomet- 
rical characterization, which requires, in principle, to consider wave propagation in different 
directions, as well as a deeper full analysis of waves. Note that here the tangle itself is not 
considered as a dynamical quantity, because it is not modified by the second sound. For this 
reason, evolution equations for the tangle are not needed here. 

The plan of this paper is the following: Section 2 is concerned with the model for helium II, in 
which the use of a pressure tensor associated to the vorticity has been considered; in Section 
3 and in Section 4 we study wave propagation in rotating frame and in presence of thermal 
counterflow respectively, pure rotation is analyzed in the general case in which a component of 
the mutual friction force parallel to the rotation axis is present; finally, in Section 5 we study 
wave propagation in simultaneous rotation and counterflow, analyzing two different situations 
about the relative direction of wave propagation with respect to the rotation vector. 

2 Evolution equations 

Many observations have shown that both thermal conductivity Ai and the relaxation time 
of the heat flux t\ in helium II are very high. As observed in \7\ their ratio — := £ < oo 
determines the velocity of the second sound, which is a heat wave propagating in the superfluid. 
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As a consequence, it is natural to use a thermodynamical theory where the heat flux q appears 
as a further fundamental field. In this way, a linear macroscopic one-fluid model of liquid 
helium II, based on Extended Thermodynamics [HE], has been formulated [7]. This model 
is able to describe the laminar flow of the superfluid both in the presence and in absence 
of dissipative phenomena and to predict the propagation of the two sounds in bulk liquid 
helium II and of the fourth sound in liquid helium flowing in a porous medium [7J, [10j-|13j, 
in agreement with microscopic and experimental data. 

In order to describe the presence of vortices in rotating helium II, in superfluid turbulence 
or in combined rotation and thermal counterflow, the use of a further additional pressure tensor 
P u , associated to the vorticity, is necessary. The simplified situations of a rotating frame and 
of pure thermal counterflow have been considered in |14j , where a constitutive relation for P w 
and its influence on the dynamics of the heat flux has been studied. 

In this work the more complex situation involving thermal counterflow in a rotating cylin- 
der, which is receiving much attention recently [15]- [21], is considered too. We start from 
a linear macroscopic one-fluid model of liquid helium II, whose fundamental fields are the 
density p, the barycentric velocity v, which is related to the two velocities of the two-fluid 
model by the relation pv = p s v s + p n v n , the temperature T and the heat flux q, related to 
the counterflow velocity V by the relation q = p s TsV (where s is the entropy of the Helium 
II). In the two-fluid model the natural variables are v s and v n , but in the experiments it is 
v and q which are directly measured. Therefore, the use of v and q appears suitable for our 
analysis. Neglecting the bulk and shear viscosity and under the hypothesis of small thermal 
dilatation (which in helium II are indeed very small), the linearized system of field equations 
for liquid helium II, in a non inertial frame, in absence of external force, is |14j : 

P^ + fe + J? + 2p(nAv) 4 = 

dT,_±_dq L= Q \ Z - 1 ) 

dt pcy dxj 

f +Cg + 2(nAq) i = (* u ) i = -(P w .q) i . 

In this system, i + 2p (ft A v) is the inertial force, £ is a positive coefficient linked to the 
second sound velocity, and: 

p = p E (p,T) and c v = (^^) (2.2) 

are the thermostatic pressure and the specific heat respectively (e is the specific internal 
energy). The effect of vortices is described by incorporating the source term P w ■ q to the 
evolution equation of the heat flux. As we will see, the expression of P w will assume different 
expressions in the different situations considered. 

Now, a small comparison between the one-fluid model and the two-fluid model could be 
useful. With the corresponding transformations between the natural variables in the one fluid 
model, v and q, and those in the two-fluid model, v s and v n , the evolution equations (12. ip of 
the one-fluid model are equivalent to those of the two-fluid model in the linear approximation 
|14j . A formal difference is found in the form of the production term in the evolution equation 
for the heat flux (12. Id ). When specified to pure rotation, this production term, as given by 
(j3.2fl . has the usual Hall-Vinen form, whereas when specified to counterflow, the production 
term, as given by (|4.ip , yields the well-known Gorter-Mellinck form. These two situations have 
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been well explored in the context of the two-fluid and one-fluid frameworks. In the combined 
situation with simultaneous rotation and counterflow, the general form of the production term 
in (12. Id ) is especially useful, as expressed in (|5.ip and f)5.2 j> . because it allows one to write 
in an explicit and appealing way the anisotropy of the tangle, whose influence on the second 
sound is one of our main concerns. Given the same geometrical conditions for the tangle - 
which here are given a priori, and whose form is probed by means of second sound — , the 
evolution equations of the two-fluid model would coincide with those of the one-fluid model. 
Thus, the dispersion relations obtained here should be valid also in the context of the two-fluid 
model. 

The one-fluid and the two-fluid models are not identical to each other. However, their 
mutual differences arise in contexts which are not relevant in the analysis presented here. 
For instance, one difference arises in the fourth sound in helium through porous media, in 
which some experimental results seem to support the one-fluid model [22]. Anyway, the two- 
fluid model could also cope with that situation provided the assumption that the superfluid 
component carries no entropy is slightly relaxed by assuming that it may carry a small but 
nonvanishing entropy. Other differences arise concerning the interaction between second sound 
and the vortex tangle. Here, we have assumed that second sound does not modify the vortex 
line density nor the geometrical structure of the tangle. If it is assumed that it may modify 
the vortex tangle, more general evolution equations would be needed, as for instance an evo- 
lution equation for the vortex line density L coupled with the rotation and counterflow, which 
have already been explored in the literature [23] . For instance, the evolution equation for L 
could be different — a generalized form of Vinen's equation with the mentioned couplings has 
been proposed and studied [16] — but this is not relevant here because an equation for L is 
not necessary in this paper, as L is taken as fixed, and its value must be found from wave 
experiments. Some other differences may appear, concerning, for instance, the possibility of 
vortex density waves at high frequencies in the one-fluid model [23] that do not arise in the 
Hall-Vinen-Bekarevich-Khalatnikov model [25]. Since in this work we are focusing our atten- 
tion to a situation in which the interaction between the second sound and the tangle does not 
distort the vortex lines nor the vortex density, the dispersion relations obtained in this paper 
by using the production terms P w (q, ft) ■ q would be also valid in the two-fluid context by 
using a production term of the form P w (v n — v s , ft) ■ (v n — v s ) in an evolution equation for 
the relative velocity V = v ra — v s . 

3 Wave propagation in rotating frame 

We generalize here the results of [14] to the case in which a small interaction between second 
sound and vortex line parallel to the rotation axis is present. In [26 1, Hall and Vinen described 
experiments of liquid Helium II in a rotating frame, showing the main effects on the propagation 
and attenuation of the second sound as a consequence of the interaction between quasi-particles 
and vortex lines: these interactions are mainly present in the planes orthogonal to the rotation 
axis. As consequence of these experiments, in [14], Jou, Lebon and Mongiovi proposed an 
expression for the production term in (12. Id ), which takes into account dissipative and non 
dissipative contributions of the interaction between quasi-particles and vortex lines, but they 
did not consider interactions parallel to the rotation axis. 

In another experiment [27], Snyder studied the component of mutual friction along the 
rotational axis, and his result, in agreement with [28], shows that this friction component is 
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very small compared with the orthogonal components but not exactly zero. In this section, 
we consider the most general case in which the axial component is included. In order to do 
that, the following vorticity tensor P w is used [T7] : 



1 



-kL r 



(b - b") ( u - nn) + b'w • Ci + 2B"£in 



(3.1) 



where U is the unit matrix, W the Ricci tensor, and B and B' are the Hall-Vinen coefficients 
[26j describing the orthogonal dissipative and non dissipative contributions while B" is the 
friction coefficient along the rotational axis. Using the Eq. (|3.ip . the production term in 
(|2.1fl ) can be expressed as [U Q2] : 



2 H 



B - B")U A (Cl A q) + B'Cl A q - 2B"Cl£l ■ q 



(3.2) 



The interest to consider spatial distribution of vortices and anisotropy of mutual friction in 
rotating container has led Mathieu et al. in [25] to analyze a more general case in which a 
parallelepipedic cavity filled of helium II rotates around an axis tilted an angle 6 with respect 
to its wall. In the following Subsection we will show that the results of the latter experiments 
can be easily explained using the general expression (|3,2p . 

Substituting the expression (|3.2p into the system (|2.ip and choosing ft = ($7,0,0), the system 
assumes the following form: 

(3.3) 

\kL r [(B - B") (- gi + qi 5 a ) - 2B" qi S a ], 
where ejyj is the Ricci tensor. 

It is easily observed that a stationary solution of this system is: 

p = Po , v = 0, T = T , q = 0. (3.4) 

In order to study the propagation of plane harmonic waves of small amplitude |30j , we put 
T = (p,Vi,T,qi), and we look for solutions of the linearized system of field equations (|2.ip of 
the form: 

r = r + fe i ^^-^, (3.5) 

where Tq = (pq,0,Tq,0) denotes the unperturbed state, f = (^p,Vi,f ', g^J small amplitudes 
whose products can be neglected, K = k r + ik s is the wavenumber, uj = uj r + iuj s the frequency 
and n = (rij) the unit vector orthogonal to the wave front. Along this paper we will assume 
that the propagating waves do not affect the vortex tangle, i.e. that they do not contribute to 
the production nor the destruction of vortices. In other terms, the waves are used to explore 
a given vortex tangle, without modifying it. If the wave amplitude is high enough, it could 
yield new contributions to the tangle. 

In the following we assume that is small, so that the term io in (|2.1b ) can be neglected. 
For the sake of simplicity, the subscript 0, which denotes quantities referring to the unperturbed 
state Tq, will be dropped out. 



8T , _J_dqj_ = „ 

dt ' pcy dxj 
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3.1 First case: n parallel to ft 

In this subsection we analyze the case in which the unit vector n orthogonal to the wave front 
is parallel to the axis of rotation, i.e. n = (1, 0, 0). Substituting (I3.5P into the linearized system 
(|3.3p and letting tj = (0, 0, 1) and t2 = (0, 1, 0) as unit vectors tangent to the wave front, the 
following homogeneous algebraic linear system for the small amplitudes is obtained: 



up + pKv\ 

OJVl 

K ' 

Ql 



-uf+ K * 











-u 



iB"KL R )qi + (KT = 



-UJV3 — 2iQv2 = 
-uiV2 + 2i£lv3 = 



(3.6) 



\kL r (B - B")) q 3 -i(2n- \kL r B') q 2 = 
\kL r {B - B")) q 2 ' 1 



+ * (20 - \kL r B') q 3 = 0. 



From the above system, it follows that longitudinal and transversal modes evolve indepen- 
dently. The study of the longitudinal modes furnishes the existence of two waves: the first is 
known as first sound or pressure wave in which density and velocity vibrate, and the second 
is known as second sound or temperature wave in which temperature and heat flux vibrate. 
Therefore, as observed in [27], when the wave is propagated parallel to the rotation axis, the 
longitudinal modes are influenced by the rotation only through the axial component of the 
mutual friction [B" coefficient). In fact, the first two equations of the system (|3.6p . give for 
the first sound V\ := ^ = Jv~p, whereas the third and fourth equation, with the assumption 
K = k r + ik s and oj real, give second sound waves with the following velocity and attenuation: 



W2 := F 



AV 2 2 k? + B^KLl *' 



W2B" ' kL r 



J R 



2Vi 



(3.7) 



where V 2 2 := is the velocity of the second sound in the absence of vortices. Therefore, the 
following fields vibrate respectively: 



cji 2 = ±kVi 



p = ip 



T 



p Y 





?i = 



^3,4 



4V 2 2 k^+B" 2 KL' 2 R 



p = 

171 = 

f = T ^ 

qi = ±pc v T 



4y 2 2 fc2+B" 2 K L| W 



On the contrary, the transversal modes are influenced by the rotation. In fact, by consid- 
ering the fifth and the sixth equation of (|3.6p they admit nontrivial solutions if and only if its 
determinant vanishes; this yields uj§§ = ±2|f2|. 
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Now, we consider the equations seven and eight of the system f)3.(5f) and, as above, we find 
the following dispersion relation: 



2Q - ^nL R B'j 



-co - -kL r (B - B") ) = 0, 



whose solutions are 



co 7fi = ±(2H - ^kL r B') - 1 -kLr(B - B") 



(3. 



(3.9) 



These transversal modes are influenced from both dissipative and nondissipative contributions 
B, B' and B" in the interaction between quasi-particles and vortex lines. 

3.2 Second case: n orthogonal to f2 

In this subsection we assume that the direction of propagation of the waves is orthogonal to 
the rotation axis, i.e. for example, n = (0, 1, 0). The unit vectors tangent to the wave front are 
ti = (1,0,0) and t2 = (0,0, 1). Under these assumptions, substituting (|3.5p into the linearized 
system (|3.3p . the following system is obtained: 

-cop + pKv2 = 



'Pp.. 
p 1 



-ujv2 + K r fp + 2ittv 3 



-COV3 



-Ulf + -^<7 2 
pcy ^ z 







-co - \kLr(B - B")) q 2 + £KT + i(2U - ±KL R B')q 3 = 
-co - \kLr(B - B")) q 3 -i(2n- \kLrB') q 2 = 



(3.10) 



— cov\ = 

(-co - iB"KLR)q\ = 0. 

In this case, the longitudinal and transversal modes do not evolve independently. The first 
sound is coupled with one of the two transversal modes in which velocity vibrates; while the 
second sound is coupled with a transversal mode in which heat flux vibrates. 

Studying the first three equations of the system (|3.10|) . we obtain a dispersion relation 
whose solutions are: 



Summarizing: 



coi = 0, 



W 2 , 3 = ±Vi 



1-4- 



9. 2 



CO. 



2,3 



(3.11) 
(3.12) 



coi = 


co 2 ,3 ~ ±KV! + o(n 2 ) 


p = tp 

V 2 = 


^ = -w * 
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The second three equations admit non trivial solutions if and only if their determinant vanishes. 
Neglecting the second-order terms in Q, the dispersion relation becomes: 



-uj — -kLr(B — B" 



-uj 



-u - ^kL r {B - B") ) - K 2 V 2 2 







(3.13) 



For uj G K and K = k r + ik s complex, one gets the solution uj/± = 0, which represents a 
stationary mode; and two solutions which furnish the following phase velocity and attenuation 
coefficient of the temperature wave: 

,2 o 



UJ 



w 2 :-- 



V4 



(B - B")kL r w 2 



W 2 2 



(3.14) 
(3.15) 



The approximated solutions to second order in 



(B-B")kL r 



w 2 ~ V 2 1 



(g - B") 2 k 2 L\ 



+ o 



are: 

{B - B"f^L R 



^ Cg - J3>L fl G f (5 - 5") 3 k 3 L| 



W 2 



UJ' 



(3.16) 
(3.17) 



Summarizing, when the direction of propagation of the waves is orthogonal to the rota- 
tion axis, the temperature wave experiences a strong attenuation, which grows with Q. The 
corresponding modes are: 



UJ4 = 



T -. 

92 = 

93 = V> 



i(2n-\nL R B') 



^5,6 



i j T/- f h (B-B") 2 k 2 L 2 r 



TnC (B-B") 2 k 2 L 2 \ , 

i(2n-fKL flJ B')T c( 1- 



93 



-," „2 ^2 



v 2 



±k r V 2 \l~'- > 



»l2 K 2 r ,2 

a l-i(B-B») K L fl 



We note that in the mode W4 = 0, only the transversal component of the heat flux is involved. 
For uj = uj r + iuj s complex and K £ K, the solutions of dispersion relation ()3.13j) are: 



0J4 
^5,6 



■-(B-B")kL r , 



±\Ik 2 v 2 



(B - B") 2 k 2 L 2 



16 



.(B - B")kL r 



R 



The first mode, with W4 = — |(J3 — B")kLr, corresponds to an extremely slow relaxation 
phenomenon involving the temperature wave and the transversal component of the heat flux: 



W4 



-UB-B")kLr 



T = 

92 = 

93 = 1p 



i(2<a-\nL R B') 
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which when f2 — > 0, converges to a stationary mode. The attenuation in uj^^q (corresponding 
to qi and ^2) is physically reasonable in view of (13, 2j) . where it is seen that for q parallel to 
fi the only component of the friction force is the axial one (related to the coefficient B"), 
whereas for q orthogonal to the vortex line (i.e. to fi) there is an attenuation dependent on 
the dissipative coefficient (B — B"). 



4 Wave propagation in presence of thermal counterflow 

In this section, we study wave propagation in presence of pure thermal counterflow in liquid 
Helium II to compare the results with those of [13] and those of Section 5. Let us consider a 
flow channel that connects two He II reservoirs (as shown in fig. 1). When a steady heat is 
applied to one end of the channel, there exists a temperature difference AT between the two 
ends. From the microscopic point of view using the two-fluid model, since only the normal fluid 
component carries entropy and heat flow, it will move away from the heat source (left reservoir) 
to the right reservoir and then give up the heat. At the same time, the superfluid component 
must counter- flow from right to left to conserve the mass. When it arrives at the left reservoir, 
part of the superfluid component will be converted to normal fluid by absorbing heat. Thus, 
a relative counterflow between the normal fluid and superfluid components is established, and 
this internal convection process is termed thermal counterflow, which is associated to the heat 
flux q through the relation q = p s TsV. 

In this case, assuming that the vortex tangle caused by the counterflow is isotropic, the vorticity 
tensor P w , as indicated in [14] . takes the following form: 



-kBLV 
3 



H 



-kBLci, 



(4.1) 



where L = 7 </ . Under this assumption, the linearized set of field equations read as: 



(dp 



di'j 



at + fe, — u 



d.r , 



n dvj I dp _ 

P at "t dxi ~ 

9T 1 _X_ dq 3 







rt pcy dxj 



k at ^ ax. 



■ 

■^KBLqi 

A stationary solution of the system (|4.2|) is [14] : 

kBL 



(4.2) 



P = Po, v = 0, T = T(x) 



3C 



-q x, q = qo 



(4.3) 



where x is the direction of the heat flux q = qo . In order to study the propagation of harmonic 
plane waves in the channel, we look for solutions of the system (j4.2j) of the form: 



p _ y 4. f e i(Kn.jXj-ui) 



(4.4) 



where Tq = (po> 0, T(x), qo), and the following homogeneous algebraic linear system for the 
small amplitudes is obtained: 



[CO 



Ujp + pKVjTlj = 

pujdi + PpKpm 
+ -^qjUj = 
+ \ikBL) qi - 





= 
= 

(Kfm 



(4.5) 



9 



The longitudinal modes are obtained projecting the vectorial equations for the small ampli- 
tudes of velocity and heat flux on the direction orthogonal to the wave front. It is observed 
that the first sound is not influenced by the thermal counterflow, while the velocity and the 
attenuation of the second sound are influenced by the presence of the vortex tangle. The 
results are: 

with p p standing for dp/ dp and: 

"•■2 - y-2\i[\ + '—^ ! =► w 2 ~v 2 (i-k 2} - 



A'. s 




'2^ 



These results generalize those of [13] where the terms in kg have been neglected. 

The transversal modes are obtained projecting the vectorial equations for the small ampli- 
tudes of velocity and heat flux on the wave front, obtaining: 

(;^ = N (4.6) 

[ (-w + %kBL) q* = 

where tt denotes the tangential plane to the wave front. The solutions of this equation are: 
uj\ = and ui2 = ^kBL. The first mode (u\ = 0) is a stationary mode. 

5 Wave propagation with simultaneous rotation and counter- 
flow 

The combined situation of rotation and heat flux (as shown in fig. 2), is a relatively new area 
of research |16j-|21|. The first motivation of this great interest is that from the experimental 
observations one deduces that the two effects are not merely additive; in particular, for q or 
Q, high, the measured values of L are always less than Lh + Lr. 

Under the simultaneous influence of thermal counterflow V and rotation speed f2, rotation 
produces an ordered array of vortex lines parallel to rotation axis, whereas counterflow velocity 
causes a disordered tangle. In this way the total vortex line is given by the superposition of 
both contributions so that the vortex tangle is anisotropic [17], [18]. Therefore, assuming that 
the rotation is along the x direction fi = (D,, 0, 0) and isotropy in the transversal (y — z) plane, 
for the vorticity tensor P^, in combined situation of counterflow and rotation, the following 
explicit expression is taken: 



'hL\1{1-D)U + D 



B"\ / \ B' B" 

1 - it ) ( u - no) + — w • Ci + 2— fiCi 

B I \ J B B 



(5.1) 



where 7 is linked to the coefficient B through the relation 7 = B/2 and D is a parameter 
between and 1 related to the anisotropy of vortex lines, describing the relative weight of the 
array of vortex lines parallel to fi and the disordered tangle of counterflow (when D = we 
recover an isotropic tangle - Eq. (|4.ip . whereas when D = 1 the ordered array - Eq. (|3.ip ). 
Assuming b = |(1 — D) + and c = ^jp-, the vorticity tensor (|5.ip can be written as: 



2b 





^ 







1-6 





1+ 








1-bJ 






7/tL { 1-6 + c V. (5.2) 



10 



Note that the isotropy in the y — z plane may only be assumed when both fi and V are directed 
along the x axis. A more general situations is yet an open topic. 

Substituting the expression (|5.2j) into the linearized set of field equations (|2.ip , it assumes the 
following form: 



dp , dv, 

+ PdXj 

dp 







dt ' pCy dXj 



jkL {2bqi5u + [(1 - b) q 2 + cq 3 ] 5 2i + [(1 -b)q 3 - cq 2 ] S 3i } 



(5.3) 



A stationary solution of this system is: 

P = Pa, v = 
T = T( Xi ) = T 



0, q = qo 

c 



(90,0,0) 
bq^SuXi. 



In order to study the propagation of harmonic plane waves, we look for solutions of (|5.3|) 
of the following form: 

r = T + fe iiKn ^- u ' t \ (5.4) 

where To = (po, 0, T(xi), qo) and T(xi) is a linear function of X{. 

Now, we investigate two different cases: n parallel to ft and n orthogonal to fi; the latter 
is the only case for which experimental data exist [6]. 

5.1 First case: n parallel to fl 

In this subsection we analyze the case in which the unit vector n orthogonal to the wave front is 
parallel to the direction of the rotation, i.e. n = (1, 0, 0). Letting ti = (0, 1, 0) and t 2 = (0, 0, 1) 
as unit vectors tangent to the wave front, the system (I5.3h for the small amplitudes (15.41) is: 

-up + Kpv\ = 



-uvx + K^p 



= 




-to 



2ij K Lb] qi + (KT = 



-lov 2 + 2iQv 3 = 
-lov 3 — 2iQv 2 = 



(5.5) 



-u 



-U) 



i^kL (1 — b)]q 2 + (2iQ — i^kLc) <?3 = 
zjkL (1 — b)] q 3 — (2iQ — ijnLc) q 2 = 

In this case the longitudinal and transversal modes evolve independently. In particular, we 
can observe that the first sound, given by the study of the first two equations of the system 
(15. 5j) . is not influenced by the presence of the vortex tangle: 

±k r Vi 



Ul, 2 
P = 



1> 
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whereas the second sound suffers extra attenuation due to the vortex tangle. The third and 
fourth equation of the system (|5.5j) admit non trivial solutions if and only if their determinant 
vanishes, obtaining in this way the following dispersion relation: 



u? + 2i-ynLbu - K 2 V 2 2 = 0. 



(5.6) 



Supposing that co is real and K = k r + ik s is complex, the dispersion relation admits the 
solutions: 



2 M t/2 



1 + a/1 + 4^L^ 



~fKLbW2 



V 2 

1 2 



(5.7) 
(5.8) 



When ft = and 6 = 1/3 the results of the Section 4 are obtained. The approximate solutions 
to second order in inLb are: 



W2 ~ v 2 1 - ^—^-a — ] + O 1 



2u< 



to 



(5.9) 
(5.10) 



Now, we study the transversal modes. The second subsystem (fifth and sixth equation) of 
the system (|5.5j) admits nontrivial solutions if and only if its determinant vanishes; this yields: 



to 2 - 4ft 2 = 0. 

The solutions of this equation are u^q = ±2101 . The respective modes are: 



(5.11) 



^5,6 


= ±2|ft| 


V3 = 




V2 = 





and they correspond to extremely slow phenomena, which, when ft — > 0, tend to stationary 
modes. 

Finally, we consider the last subsystem (equations seven and eight), whose dispersion 
relation is: 



u} 2 + 2ijKL(l -b)uo+ - (jkL(1 - b)) 2 - 4ft 2 + 4 7 ftftLc - (jkLc) 2 
which admits the following exact solutions: 

ujj t s = =t (2ft — 'jkLc) — ijuL (1 — b) . 
The corresponding modes are: 

ojj ; s = =t (2ft — 7«Lc) — i^kL (1 — b) 







(5.12) 



(5.13) 



93 = ip 
q 2 = ±iip 
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From (15. 71) . (I5.8P and (I5.13P one may obtain the following quantities L, b and c: 



L 



-U S W2 + V^fc, 



■uj r W2 + 2f2u>2 



(5.14) 



where we have put uj = u r + iu s . 

The results of this section, from the physical point of view, imply that measurement in a single 
direction are enough to give information on all the variables describing the vortex tangle. 

5.2 Second case: n orthogonal to fi 

Now we assume that the direction of propagation of the waves is orthogonal to the rotation 
axis, i.e. for example, n = (0, 1, 0). The unit vectors tangent to the wave front are ti = (1, 0, 0) 
and t2 = (0,0, 1). Under these assumptions, the homogeneous algebraic linear system for the 
small amplitudes is: 



-up + Kpv2 = 



T> 

2iQ,V2 



-uv 2 + K^jp + 2iriv 3 = 



-uv 3 



-uq 2 + QKT + 2i£lq 3 = I^kL [(1 - b)q 2 + cq 3 ] 
-uq 3 - 2iQq 2 = ijkL [(1 - b)q 3 - cq 2 ] 



(5.15) 



-uv\ = 

— u — i"fK2Lb] qi = 

In this case the longitudinal and the transversal modes not evolve independently. In particular, 
the first sound is coupled with one of the two transversal modes in which velocity vibrates, 
while the second sound is coupled with a transversal mode in which heat flux vibrates. 
As in the previous subsection, the first subsystem (first three equations) of the system (|5.15p . 
admits non trivial solutions if and only if its determinant vanishes: 



u u 



[u 2 - AU 2 - K 2 p p ] = 0. 



(5.16) 



The solutions of this equation, as also the corresponding modes, are the same to the case of 
pure rotation (see equations (I3.11j) - (l3.12p ). 

The second subsystem (fourth and fifth equations), has the dispersion relation: 



( — u - hkL(1 - 6)) [u (-u - i~fKL{l - b)) + K 2 V 2 2 ] + u (2«Q - i-fKLc) 2 = 

Assuming u S 3ft and K = k r + ik s , one obtains the following two equations: 

-u 3 + 7 VL 2 (1 - b) 2 u + 4n 2 u + -f 2 K 2 L 2 c 2 u - AjkLcQu + 
+k 2 .V 2 2 u - k 2 V 2 2 u - 2 7 kL(1 - b)k r k s V 2 2 = 0, 



0. 



(5.17) 



-27«L(1 - b)u 2 + 2k r k s V 2 2 u + 7«L(1 - b)(k 2 - k 2 )V 2 2 



0. 



(5.18) 
(5.19) 
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In the hypothesis of small dissipation (k 2 kg), from (I5.19P one obtains: 



k s = jkL(1 — b) 



which substituting in (|5.18p . yields: 



2w\ 



v 2 



2w 2 V 2 2 



(5.20) 



(2fi - -ikLcY - 7 Vi/(i - by u 1 - KVjuj* - 7 Vz/(i - byvfK = o. (5.21 



Putting A = — 
account that w 2 



(2fl- 1K Lc) 2 - 7 VL 2 (l-6) 2 
= f-, the Eq. (ET211 becomes: 



and B = — ^ 2 k L (1 — by and taking into 



whose solutions are: 



1 + 



A 



oj- 



V4 



B 



(5.22) 



Wo 



0, and 



u?-B 



u 2 + A 



V? 



1 



(2Q-jKLcy 
w 2 +7 2 k 2 L 2 (1-6) 2 



(5.23) 



We can remark that the coefficients A and B are negative and that w\ > V 2 2 because to 2 + A < 
uj 2 — B and, in particular, w 2 = V 2 for f2 = 1K Lc . Now, studying the transversal modes, i.e. 
the third subsystem (equations seventh and eighth), we obtain = 0, which corresponds to 
a stationary mode, and: 

uj s = -i^K2Lb. (5.24) 

Summarizing, also in this case measurement in a single direction are enough to given infor- 
mation on all the variables describing the vortex tangle, namely L, b and c, from equations 
(15^01) . (15^31) and (IOTP : 



L 



4k s w 2 V 2 2 - u s (2w 2 - Vo 2 ) 



(2w 2 



{2w 2 



V 2 ) 7 K 



4k s w 2 V 2 2 

m(2wl 



(2w 2 - V 2 ) ' 



(5.25) 



V 2 ) - V(l - V 2 )(4k 2 (2w 2 - V 2 ) 2 + W^Vj) 



4k s w 2 V 2 2 



U! s (2w 



V?) 



where we have put uj% = iuo s and u s = 2 r ynLB. 

In this subsection we have analyzed wave propagation in the combined situation of rotation 
and counterflow with the direction n orthogonal to Q. In [6] Swanson et al. experimented 
the same situation, but they didn't represent the attenuation neither the speed of the second 
sound but only the vortex line density L as function of and V. Therefore, it is unknown how 
they plotted these graphics, which the hypothesis were made and what was the anisotropy 
considered. Instead, the results of these two subsections allow to know the spatial distribution 
of the vortex tangle simply by performing experiments on waves propagating orthogonally to 
Q (equations (|5.14p ) or parallelly to Q (equations (|5.25p ). 
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6 Conclusions 



In this work we have studied the propagation of waves (longitudinal density and temperature 
waves, longitudinal and transversal velocity and heat waves) in turbulent superfluid helium in 
the three situations: rotating frame, thermal counterflow, and simultaneous thermal counter- 
flow and rotation. 

From the physical point of view it is interesting to note that our detailed analysis in Section 5 
shows that, in contrast to which one could intuitively expect, measurements in a single direc- 
tion are enough to give information on all the variables describing the vortex tangle, namely 
L, b and c, for instance, from one of (j57Tj) - (jOj) and (pU3]) or of (l5T2T)j) - (|5T23"1 ) and (JEM]). This 
is not an immediate intuitive result. Future analyses of work along this direction could be, for 
instance, to consider that and V have arbitrary directions, i.e. that they are not parallel to 
each other, in which case (|5.2p would not be sufficient to describe the vortex tangle, because 
no isotropy in the y — z plane could be assumed. 

Another topic could be to assume that the external waves produce vibrations in the vortex 
lines, without creating nor destroying them. An example of that is the work of Barenghi 
et al. [28J . A more general possibility would be to consider that nonlinear effects of the 
external waves create and destroy new vortices. Yet another topic would be to consider what 
happens with waves whose wavevector A become short enough to be comparable with the 
average vortex separation, of the order L _1//2 . In this case, one could study nonlocal effects in 
the vortex |33|. [34] . The first mentioned application could be carried out within the existing 
physical model, at the expenses of more cumbersome calculations. In contrast, the other three 
applications need more progress in the basic physical understanding of the problem. 

Acknowledgments 

We acknowledge Prof. D. Jou (Departament de Fisica, Universitat Autonoma de Barcelona) 
and Prof. M.S. Mongiovi (Dipartimento di Metodi e Modelli Matematici, Universita di 
Palermo) for the enlightening discussions useful to the deepening of these arguments. This 
work is supported by MIUR of Italy under "Azioni Integrate Italia-Spagna, anno 2005". R.A. 
P. is supported by an " Assegno di ricerca delklstituto Nazionale di Alta Matematica 'F. Severi' 
(INdAM)" of Italy and M.S. is supported by an "Assegno di ricerca MIUR" of Italy. 

References 

[1] L. Tisza, Nature 141 (1938) 913. 

[2] L.D. Landau, J.Phys. 5 (1941) 71. 

[3] J.T. Tough, Prog. Low Temp. Phys. 8 (1982) 133. 

[4] R.J. Donnelly, Quantized Vortices in Helium II, Cambridge University Press, Cambridge, 
UK, 1991. 

[5] R.T. Wang, C.E. Swanson, R.J. Donnelly, Phys.Rev.B 36 (1987) 5240. 

[6] C.E Swanson, C.F. Barenghi, R.J. Donnelly, Phys.Rev.Lett. 50 (1983) 190. 

[7] M.S. Mongiovi, Phys. Rev. B 48 (1993) 6276. 



15 



[8] D. Jou, J. Casas- Vazquez, G. Lebon, Extended Irreversible Thermodynamics, Springer- 
Verlag, Berlin, 2001. 

[9] I. Miillcr, T. Ruggeri, Rational Extended Thermodynamics, Springer- Verlag Berlin, 1998. 
[10] M.S. Mongiovi, R.A. Peruzza, Zangew.Math.Phys. 54 (2003) 566. 
[11] M.S. Mongiovi, R.A. Peruzza, I. J. Nonlinear Mech. 39 (2004) 1005. 
[12] M.S. Mongiovi, R.A. Peruzza, Math. Comp. Model. 38 (2003) 409. 

[13] M.S. Mongiovi, R.A. Peruzza, D. Jou, Recent research developments in Physics, 
Transworld Research Network, (2004), 1033. 

[14] D. Jou, G. Lebon, M.S. Mongiovi, Phys. Rev. B 66 (2002) 224509. 

[15] C. F. Barenghi, R.J. Donnelly, W.F. Vinen, Quantized Vortex Dynamics and Superfluid 
Turbulence, Springer-Berlin, 2001. 

[16] D. Jou, M.S. Mongiovi, Phys. Rev. B 69 (2004) 094513. 

[17] D. Jou, M.S. Mongiovi, Phys. Rev. B 72 (2005) 144517. 

[18] D. Jou, M.S. Mongiovi, Phys. Rev. B 74 (2006) 054509. 

[19] M.S. Mongiovi, D. Jou, Phys. Rev. B 72 (2005) 104515. 

[20] M. Tsubota, C.F. Barenghi, T. Araki, A. Mitani, Phys. Rev. B 69 (2004) 134515. 
[21] M. Tsubota, T. Araki, C.F. Barenghi, Jour. Low Temp. Phys. 134 (2004) 471. 
[22] M.S. Mongiovi, Physica A, 291, (2001) 518. 
[23] M.S. Mongiovi D. Jou, Phys. Rev. B, 75, (2007) 024507. 

[24] D. Jou, M.S. Mongiovi, M. Sciacca, "Vortex density waves in a hydrodynamical model of 
superfluid turbulence", Phys. Lett. A, in press 

[25] K.L. Henderson and C.F. Barenghi, Theoret. Comput. Fluid Dynam. 18 (2004) 183. 

[26] H.E. Hall, W. F. Vinen, Proc. R. Soc. London A238 (1956) 215; Proc. Roy. Soc. A238 
(1956) 204. 

[27] H.A. Snyder, Physics of Fluids 6 (1963) 755. 

[28] C.F. Barenghi, M. Tsubota, A. Mitani, T. Araki, J. Low Temp. Phys. 134 (2004) 489. 
[29] P. Mathieu, B. Placais, Y. Simon, Phys. Rev. B 29 (1984) 2489. 
[30] J.Whitham, Linear and Nonlinear Waves, New York, Wiley, 1974. 

[31] I.M. Khalatnikov, An Introduction to the Theory of Superfluidity, New York, Benjamin, 
1965. 

[32] Y.A. Sergeev, C.F. Barenghi, J. Low. Tern. Phys. 127 (2002) 203. 

[33] M.S. Mongiovi, D. Jou, Phys. Rev. B 71 (2005) 094507. 

[34] M.S. Mongiovi, D. Jou, J. Phys: Condens. Matter 17 (2005) 4423. 



16 




Fig. 1. Counterflow container configuration. 




Fig. 2. Rotating counterflow container configuration. 
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